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Abstract 

The applied method of slowly varying amplitudes of the electrical 
and magnet vector fields give us the possibility to reduce the nonlin- 
ear vector integro-differential wave equation to the amplitude vector 
nonlinear differential equations. It can be estimated different orders 
of dispersion of the linear and nonlinear susceptibility using this ap- 
proximation. The critical values of parameters to observe different 
linear and nonlinear effects are determinate. The obtained amplitude 
equation is a vector version of 3D-I-1 Nonlinear Schredinger Equation 
(VNSE) describing the evolution of slowly varying amplitudes of elec- 
trical and magnet fields in dispersive nonlinear Kerr type media. We 
show that VNSE admit exact vortex solutions with classical orbital 
momentum £ = 1 and finite energy. Dispersion region and medium 
parameters necessary for experimental observation of these vortices, 
are determinate. 

PACS 42.81.Dp;05.45.Yv;42.65.Tg 

1 Introduction 

In present day, there are no difficulties for experimentalist in laser physics 
or nonlinear optics to obtain picosecond or femtosecond optical pulses with 
equal duration in x, y and z directions. The problems with so generated light 
bullets arise in the process of their propagation in a nonlinear media with 
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dispersion. In the transparency region of a dispersive Kerr type media, as it 
was established in jH I2l E] , the scalar paraxial approximation (no dispersion 
in the direction of propagation), is in very good accordance with the exper- 
imental results. The paraxial approximation, used in the derivation of the 
scalar 2D+1 nonlinear Schredinger equation (NSE), do not include a small 
(in first or second order of magnitude) second derivative of the amplitude 
function in the direction of propagation. The including of this term does 
not change the main result dramatically: in the case of a linear propagation 
(pulses with small intensity), as generated optical bullets at short distance 
are transformed in optical disks, with large transverse and small lengthwise 
dimension. Some of the experimental possibilities, this small term to become 
important, were discussed in In Ref. [H] it was shown, that this second 
derivative in the direction of propagation term is in the same order and with 
the same sign as the others only in some special cases: optical pulses near 
the Langmuir frequency or near some of the electronic resonances. In these 
regions the sign of the dispersion is negative and the scalar 2D+1 NSE be- 
comes 3D+1 NSE ones jSl 13 IE]- Propagation of optical bullets under the 
dynamics of 2D+1 and SD-fl NSE are investigated also in relation to dif- 
ferent kind of nonlinearity CD]. A generation of a new kind of 2D and 
3D optical pulses, so called optical vortices, has recently become a topic of 
considerable interest. Generally, the optical vortices are such type of optical 
pulses, which admit angular dependence of electrical field or helical phase 
distribution. The electrical field or intensity is zero also in the center of the 
vortices. The original scalar theory of optical vortices was based on the well 
known 2D-I-1 NSE p!T|ll2[ [T^. In a self focusing regime of propagation can be 
generated optical rings, but they are modulationally unstable [13 HE] ■ One 
alternative way of stabilizing of optical vortices in 2D and 3D case, using 
saturable [T7| HH] or cubic- quintic [121 12D] nonlinearity, was discussed also. 
From the other hand, the experiments with optical vortices show, that the 
polarization, and the vector character of the electric field, play an important 
role in the dynamics and the stabilization of the vortices pji . To investigate 
these cases we are going to a vector version of 3D+1 NSE. In Section 2 we 
derive this vector version of 3D+1 NSE (VNSE), using also the dispersion 
in the direction of propagation and envelope approximation in the standard 
way as it was constructed in [21 1^- This vector generalization allows us to 
find exact vortex solutions with spin 1=1 and finite energy. The question 
for shape and finiteness energy of these vortices is crucial, and will discuss 
largely in Section 5. 
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2 Maxwell's equations for a source-free, dis- 
persive, nonlinear Kerr type medium 

The Maxwell's equations for a source-free dispersive media with Kerr non- 
linearity are: 



V • D = 0, (3) 
V ■ B = V ■ H = 0, (4) 

B = H, D = Pu,, + 4TTpnHn, (5) 

where E and H are the electric and magnetic intensity fields, D and B are the 
electric and magnetic induction fields, Pun, Pniin ai'c the linear and nonlinear 
polarizations of the medium respectively. The linear electric polarization of 
an isotropic, disperse medium can be represented as: 

t 

Piin^ J {5{t-r)+A7rx^^^(t-r))E{T,x,y,z)dr^ 

— oo 

t 

j £o{t-T)E{T,x,y,z)dT, (6) 

— oo 

where x^^^ are the linear electric susceptibility and the dielectric con- 

stant. In the following, we will study such polarization (linearly or circularly 
polarized light) , where the nonlinear polarization in the case of one carrying 
frequency can be expressed as: 
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—CO —CO — oo 

X (-E(ri, X, y, z) ■ E*{t2, x, y, 2:)) E{t3, x, y, z)dTidT2dT3. (7) 

It is important to point here the remark of Akhmanov at all in |22] , that the 
method slowly varying amplitudes of electrical and magnet fields applied to 
such nonstationary linear and nonlinear representation is valid as well when 
the optical pulse duration of the pulses to is greater than the characteristic 
response time of the media Tq {to » Tq), as when the time duration of the 
pulses are equal or less than the time response of the media (to < Tq). This 
possibility is discussed in the process of deriving of the amplitude equations. 
Taking the curl of equation ((T)) and using Q and (0), we then obtain: 



ae-v(v.e)-1^ = o, (8) 

where A = is the Laplasian operator. Equation (jH)) is derived without the 
use of the third Maxwell's equation. Using equation (jSJ and the expression 
for the linear and nonlinear polarizations © and (|7j), it can be estimate the 
second term in the equation (jHJ for arbitrary localized vector function of the 
electrical field. As is shown in for such type of function V ■ -E = and 
equation (8) is written: 



A^-i^^O. (9) 



c 



We define a complex presentation of the real electrical field by the relation: 



E {x, y, z,t) = A (x, y, z, t) — (exp {i{kQZ - uqI)) - c.c) , (10) 

where A' and ujq and are the real vector amplitude, the optical frequency 
and the wave number of the optical field respectively. The real vector ampli- 
tude A' is represented also by a complex vector amplitude field: 
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1' {x, y, ^,t) = ^. (^A{x, y, z, t) - A* {x, y, z, t)) , (11) 

where A is the complex ampHtude of the electrical field. In this way our real 
electrical field is presented with four complex fields of kind: 



E {x, y, z, t) = — - (^Aexp {i{koZ — uJot)) + A* exp {—i{koZ — oj^t 

+ ^ (^A exp {—i{kQZ — Uot)) + A* exp {i{kQZ — UqI))^ . (12) 

This special kind of presentation of electrical field is connected with shape 
and finiteness energy of vortex solutions, and it is discussed in Section 5. Here 
are used also the Fourier representation of the complex amplitude function 

— * 

A and their time derivatives to second order: 



A{x,y,z,t)= j A{x,y, z,uj — (jJq) exp {—i {(jJq — uj)t)duj, (13) 



dA{x,y,z,t) 
dt 



-co 

-I 

-oo 



i j {uj - ujq) A{x,y,z,uj - ujo) 

X exp (— i (a; — ujq) t)du!, (14) 



+ 00 

j {u - uof A{x,y,z,uj - uo) 

— oo 

X exp (i (uj — ujq) t)duj. (15) 

The principe of causality request the next conditions on the responce func- 
tions: 

£(i-T) =0; x^'U^-^i>^-^2,t-T-3) = 0, 

t-T < 0; i -Ti < 0; i = 1,2,3. (16) 



d^A{x,y,z,t) 
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That because we can prolonged the upper integral boundary to infinity and 
to use standard Fourier presentation [H]: 



t 



J eo{T — t) exp {iujT)dT = J £o{t ~ t) exp {iujT)dT, (17) 

-oo — oo 

t t t 

III ('^ ~ ^1' ~ ^2,t - 1"3) '^1"1'^1"2C?T3 = 

— oo — oo — oo 

+00 +00 +00 

/ / / ^^'^ {t - Tut - r2,t - Tg) dndT2dT3. (18) 



-oo — oo — oo 



The spectral presentation of linear optical susceptibility £:'o('^) is con- 
nected to the nonstationary optical response function by the next Fourier 
transform: 



80(10) exp {—iut) = / eo{t — t) exp {—iLUT)dT. (19) 



Similar is the expression for the spectral presentation of the non-stationary 
nonlinear optical susceptibility x^^^: 



hoo +00 +00 



X'^'\uj = 2lo-u) exp i-iujt) = J J J X^'^ {t - n,t - r2,t - rg) 

—00 —00 —00 

X exp {-i {uj{Ti +T2 + Ts)))dTidT2dT3. (20) 

From (jni), (frnjl . (fT!?jl (fTTjl and (fnijl for the linear polarization of one of the 
complex fields P/j„ A exp {i{kQZ — cut)) is obtained: 



t 



Plin = I ^0 (t - t) (K^Z - UJoT)) J A{x,y,Z,UJ ~ (Jo) 

—00 

X exp {—i{uj — ujQ)T)dujdT 
exp {ikoz) J Ai (x, y, z,uj — cuq) J eo{t — r) exp {—iujT)dTduj. (21) 



00 —00 



+00 



—00 —00 
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The second integral in (|2ip is simply equation (|19|). and the linear polariza- 
tion can be written: 



Pii^{x,y,z,t) = exp{i{koz - ujot)) / eo{uj) A{x,y, z,uj - ujq) 

X exp (— i {uj — ujq) t)duj. (22) 



It is important to point, that the similar expressions as ()22|) are written also 
for the polarization of the complex fields of kind: 



and 



lin 



lin 



y4exp {—i{koZ — cut)) ; P^^ A* exp {i{koZ — out)) 



Aexp {i{koZ — ujt)) ; P^^^ A* exp {—i{kQZ — ut)) 



(23) 



(24) 



The above procedures are used also to the nonlinear polarization. When the 
third harmonic term is neglected the result is: 



P''^'''{x,y,z,t) = -exp{i{koz-uJot)) j x^^\uj) exp {-t{uj - uJo)t) 

— oo 

X \A{x, y, z,uj — iJo)\'^A{x, y, z,tu — uJo)dtu. (25) 

Multiply (I22l), ^ by 1/c^ and using the second derivatives in time, we 
obtain: 



1 d^Piin (x, y, Z, t) 

.2 



— exp {i {k^z — ujQt)) 



A{x,y,z,uj - ujq) 



X exp (— i [uj — cuo) t)(ico',(26) 



— = -exp{t{koz-uJot)) J 



x\A{x, y, z,u — uJo)\'^A(x, y, z,uj — cuq) exp {—i{uj — uJo)t) dcu. (27) 
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The spectrum of the amphtude function is restricted, by writing the wave 
vectors: 



9 ' 



and 



near the carrying frequency in a Taylor series: 



kL (^) = J = ^ (^o) + - ^o) 



(28) 



k. ^ f-^^^ (29) 



+-^^(..-^0) +..., (30) 



k'nl (^) = = (^O) + ^ g'J^ (^ - ^O) + ... (31) 

It is convenient to express the nonhnear wave vector by hnear wave vector 
and nonhnear refractive index: 



j2f N Q 2-^(3)/ N/2 ^O^(^o) 37rx(^)(^o) ,2 .oo^ 

kniM = ^'^^oX^ 'M/c = — ^ — = kon2, (32) 



where 



^2^0 = — 7^—^ — . (33) 

is the nonhnear refractive index. We note that any approximation of the 
response function no used. There is only one requirement of the spectral 
presentations (fT^ and (j^UI) of the response functions: to admit first and sec- 
ond order derivatives in respect to frequency (to be smooth functions). So, 
this method is not limited from the time duration of the response func- 
tion. This give us the possibility to apply it also when the half-max of the 
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pulses are in order of the time duration of the response function (femtosec- 
ond pulses). The restriction is only in respect of the relation between the 
main frequency ujq and time duration of the envelope functions to determinate 
from the relations ()30|) and (jHT|) (conditions for slowly varying amplitudes). 
Putting eqn. (jHn|) and (jHTj) in (f^ and (jTTj) respectively and keeping in mind 
the expressions of time derivative of the amplitude functions (I14|) - (jl5p for 
the electric field, the second time derivative of the linear polarization of the 
optical component is represented in next truncated form: 



1 d'^Piin {x, y, t) 

C2 dt^ 



2iko dA 
V dt 



+ hK + 



1 \ d^A 



+ ... 



X exp (i(A;o-2 — cuoi^))- (34) 



Similar result is obtained also for the linear polarization of the terms of kind 
(\2'i\ and (j24j) . The nonlinear polarization therm is: 



47r d'^Pnlin {X, y, Z, t) 
C2 dt^ 



( 



-kln2 



A 



A + 2ik 



dk. 



nl 



d 



A 



A 



nl ' 



did 



dt 



+ 



X exp {i{kQZ — uj^t)) 



I 



n2 



-kQn2 A A + \ 2iko h ik, 



V 



' a7 



d 



A 



A 



dt 



+ ... 



X exp {i{koz — ujQt)) , (35) 



where v = l/{dkiin {ujq) /duo) and /cq are the group velocity, and dispersion. 
From the wave equation ©, using JHI), (jSH), ^ and (jlSl) , (El) , the next 
two equations of complex vector amplitude and it complex conjugate are 
obtained: 



dA dA ( kov dn2 ' 



d 



A 



V 



2kn 



V 

-AA-- 



1 \ d^A 



k^v"^ J dt 



+ 



dt 



A 



A. 



(36) 
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The equations and (jHTj) are equal and solving for A we can find 
also A* and respectively the real amplitude vector field A'. We investigate 
the case when: 



(38) 



a condition that can only be fulfilled for materials possessing negative dis- 
persion. The term in this case is neglected. Applying a Galilean trans- 
formation to the vector amplitude equation ()36|). where the new reference 
frame moves at the group velocity, t' = t; z' = z — vt (the primes are missed 
for clarify) , we obtain our final vector amplitude equation: 



3 Norming VNSE. Exact vortex solutions 



In this section we will obtain solutions in the case of amplitude vector equa- 
tions of one carrying frequency. The case of two carrying frequencies is 
discussed in jH]. The case of three and more localized optical waves on dif- 
ferent frequencies with additional conditions on the frequencies (parametric 
vortices) will be discussed in a next paper. Defining the rescaled variables: 




(39) 



where A 



± 



A = AqA'; X = rox'; y = rov'; z = rQz'; t = t^t', 



(40) 
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and constants: 



a = 2korQ/tov; f3 = v fc"fco; 7 = ^o^o^a |^o| ; 

7i = koron2 \Aof ; 72 = -koron2 \Aof , 

c 



(41) 



equation can be transformed in the following (the primes are not writ- 
ten): 



— z 



d 



OA , 



A 



'a]\ 



\ 



d'^A 

A±A-/?^+7 A 



A = 0. 



(42) 



The dispersion therm (the second derivative in 'z' direction of the amplitude 
function) in transparency region of a media is usually one or two order of 
magnitude smaller than the diffractive therm (transverse Laplasian). That 
because usually the linear (dispersion) parameter is /5 ~ 10~^. There is 
possibility to reach (3 = —1 only in some special cases, near to Lengmuire 
frequency in cold plasma, high-frequency transparency region of dielectrics 
and near to some of the electronic resonances j3] . The constant a has a value 



of a ^ 10' 



(a 



rofco) if the of slowly varying approximation is used. When 



the nonlinear constant admit typical, critical for self-focusing regime value 
7 = 1, the constants 71 = 72 ~ 10~^ are small. We point here that the effects 
of asymmetry of the pulses (asymmetry of their spectrum) due to nonlinear 
addition to the group velocity presented by second therm in equation is 
substantial when 71 = 72 ~ 1 or 7 ~ 10^. These terms depends from intensity 
of the fields and when 7 ~ 10^ or n2 |Ao|^ ~ 10~^ effects of self-steepening of 
the pulses, connected also with considerable self-focusing can be estimated. 
Such type of experiments are provided with 80 femtosecond pulses and 
intensity Jq ~ lO^'^Wt/cm'^. This correspond to the parameters discussed 
above (7 ~ 10^ or n2 {Aq^ ~ 10~^). In this paper is investigated the case, 
when /3 = — 1 (negative dispersion) and 7 = 1. Neglecting the small terms 
the equation becomes: 



Now we going to the next step - the possible polarization of the ampli- 
tude vector field Recently, this problem was discussed largely ^^ ^^ 
and was pointed the different between the polarization of the optical waves 
with spectral bandwidth (slowly varying amplitudes), and the polarization of 
monochromatic fields. In case of monochromatic and quasi-monochromatic 
fields the Stokes parameters can be constructed from transverse components 
of the wave field [21]. This lead to two component vector fields in a plane, 
transverse to direction of propagation. For electromagnetic fields with spec- 
tral bandwidth (our case) the two dimensional coherency tensor cannot be 
used and the Stokes parameters cannot be found directly. As it was shown by 
T.Carozzi and all in [2^], using high order of symmetry (SU(3)), six indepen- 
dent Stokes parameters can be found. This corresponds to three component 
vector field. Here is investigate this case. The increasing of the spectral 
bandwidth of the vector wave, increasing also the depolarization term (com- 
ponent, normal to the standard Stokes coherent polarization plane). The 
amplitude vector function of electrical field A is represented as sum of three 
orthogonal linearly polarized amplitudes: 



A{x, y,z,t) = J2 J^j ^) 6xp (i^t) . 



(44) 



In a Cartesian coordinate system, for solutions of the kind of ()44|) , the vector 
equation is reduced to a scalar system of three nonhnear wave equations: 



aQA^ + AA^+ 



A^ 



A' 



0; / = x,y,z. 



(45) 



As it was pointed in |3J, we may choose to express each of the components Ai 
in spherical coordinates Ai{r,6,ip) of the independent variables i = x,y,z. 
In this way the linear unique (polarization) vector of each of the components 
hold down. The system ()45|) in spherical variables is: 



aQA' + ArA' + -Ae,^A' + ^ 



A^ 



A^ = 0; 1 = x,y,z, 



(46) 



where 



12 



sm 9 o9 \ o9 J sm 9 o^p^ 



are the radial and the angular operator respectively and r = \/ x"^ + y"^ + 
9 = arccos ^ ip = arctan^, are the moving spherical variables of the inde- 
pendent variables x, y, z. The system of equations pUj) are solved, using the 
method of separation of the variables. We present the components of the 
field product of a radial and an angular part: 



A\r,9,cp) = Rir)Yii9,cpyj = x,y,z, (49) 
with the additional constraint on the angular parts: 

\Y^{9,^)\^ + \Yy{9,^)\^ + %{9,^)f = const. (50) 
Multiplying equation ()46|) by and bearing in mind the constraint ex- 



pressed in (lHn|) . we obtain: 

^ + {an + \Rf) = = m + 1), (51) 

where £ is a number. Thus the following equations for the radial and the 
angular parts of the wave functions are obtained: 



ArR + anR+\RfR-^^^-—^R = (52) 



i{i + 1) 



Ae,^Yi + i{i + 1)1^ = 0; / = x, y, z. (53) 

Equations ()52|) - (j3nj) shows that the nonlinear term occurs only in the radial 
components of the fields. As for the angular parts we have the usual linear 
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eigenvalue problem. The solutions of equations (f^ of the angular parts are 
well known; each of them has exact solutions of the form: 



Yi = Y[ 




2£+ (£ 



47r(£ + m)\ 



■ (cose) exp (imip), (54) 



where P/" are the Legendre's functions for a discrete series of numbers: i = 
0, 1, 2....; m = 0, ±1, ±2, .. and \m\ < i. Returning to equation it is seen 
that separation of variables for spherical functions, which satisfy condition 
flMH) . is possible only for / = 1: 



Yx = Y-^ ^ = sin 9 cos (p 
Yy = Y^ = sin 9c sin Lp 
Y, = Y^ = cos 9, 



(55) 



or another appropriate combination of axes. By choosing one of these angular 
components for each of the field components we see that the eigenfunctions 
fl^ are solutions to the angular part of the set of equations fl^Hj) . It is 
straightforward to show that the radial part of the equations fl52|) admit "de 
Broglie soliton" solutions |EZI of the form: 



R=V2 



exp (j,y/aQr 



r 



(56) 



In this way, we prove the existence of vortex solutions with angular momen- 
tum / = 1 in a Kerr type medium. The solutions of the vector amplitude 
equation in a fixed basis are: 



exp 



sin 9 cos if exp (iQt) , 



(57) 



A^ = V2' 



exp (i 



'aQr 



sin 9 sin ip exp [iQt) , 



(58) 
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A, = V2 ^ ^ cos e exp {iVlt) 



(59) 



The equation is normed and this lead to the next normed constants: 



aQ = l: Q =~ 10" 



(60) 



Now we can turn bask to the amphtude of real solutions A' , using equation 



a: 



2i 



{A^ — Al] 



V2 



sin (^\/ aflr + Qt 



sin 9 cos ip, 



(61) 



A 



2i v^'^ " 



V2 



sin (^V aVtr + Vtt 



sin^sincyj, 



(62) 



^sin ( a/ af2r + fit 
\/2 ^ 



cos 6*. 



(63) 



The condition = separates the nonlinear Maxwell's system in two wave 
equation system, one nonlinear, of the electric field, and another nonlinear 
wave equation of the magnetic field. We define again a complex presentation 
of the real magnet field by the relation: 



H (x, y, z, t) = C (x, y, z, t) — (exp {i{kQZ - uj^t)) - c.c) 

2i 



(64) 
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where C' and ujq and fco are the real vector amphtude, the optical frequency 
and the wave number of the optical field respectively. The real magnet vector 
amplitude C is represented also by a complex vector amplitude field: 

C' {x, y, z,t) = ^. (C (x, y, z, t) - C* (x, y, z, t)) , (65) 

where C is the complex amplitude of the magnet field. Applying similar 
procedures to those who were made for electric, namely; using the Fourier 
representation of the amplitude functions C and their time derivatives to 
second order, as were done for the electric, we obtain the vector equation of 
slowly varying amplitudes of the magnetic field C: 



. dC , ^ d'C 



A 



'C = 0, (66) 



if the condition: 



A 



X A = 0, 



(67) 



is satisfied. For solutions of the amplitudes of electrical field like (j^ - lj^UI) . 
the condition (jU7j) is satisfied. Equation (jMjl has a localized solution of the 
same kind (j37jl - (j^ . as the amplitude functions of the electric field, but with 
opposite phase: 



^exp ( ailr 
V2 ^ 



sin 6 cos (f exp {—iQt) 



(68) 



exp 



sin 6 sin exp {—iQt) 



(69) 



^exp ( iv aQr] 

C, = V2 ^- ^ cos 9 exp {-iQt) . (70) 

The real dynamics of localized fields can be understood only by investi- 
gating both equation (jl^ and (jMjl .The solutions of electrical (jFTjl - and 
magnet (jUHjl - (fTOj) fields are with opposite phase, and that because the cor- 
responding electromagnetic localized wave oscillating in coordinate system 
moving with group velocity. 
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4 Experimental conditions 



The results presented above require the conditions relating the linear param- 
eter (5 and the nonlinear parameter 7 to satisfy: namely (3 = —1 and 7 = 1. 
The relation 7 = 1 is a typical critical value for starting of self-focusing 
regime of a Gausian pulse. The nonlinear parameter written again is: 

7 = ro/con2|v4o|^ = a^nsl^ol^ = 1. (71) 
The constant a = r^ki in the optical region ranges in value from: 

a = 10^ - 10^ (72) 

Using this range in the condition (|71|). we obtain a required nonlinear refrac- 
tive index change: 

nal^ol^ = 10"^- 10^^ (73) 

It was shown, using relations (j4ip . that if we investigate near to the criti- 
cal value of self-focusing 7 = 1, then 71 ~ 72 10~^ — 10~^. This is one 
important result: The nonlinear addition to the group- velocity therm in am- 
phtude equation is relatively small, even for (femtosecond) pulses with 
time duration equal or little to characteristic time response of the nonlinear 
media. 

Using the condition of linear parameter 13: 

f3 = kov^k" = -1, (74) 

we see that the vortex propagation take place only in the negative-dispersion 
region. In [H] we find, that the constraint, given in equation (|74p correspond 
to the next two experimental situations: 

1. Cold plasma near the Langmuir frequency. 

2. A region near an electronic resonance in an isotropic medium. 

There is also a equivalence in high-frequency region between the spec- 
tral presentation of dielectric susceptibility of a cold plasma and this of a 
dielectric media. The expression is equal wish precise a (dielectric) constant. 
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Near these frequencies the dispersion parameter increases rapidly and admits 
values: 

I A;" |~ IQ-'^^ - iQ-'^^i^ This lead to the fact, that the dispersion term 

I I cm ' -t^ 

(the second derivative in z direction of the amplitude function) in normed 
amplitude equations of electric (jl^ and magnet field is of the same 
order of the diffractive term (transverse Laplasian). 



5 Finiteness of the energy for the vortex so- 
lutions 

Now we come to the crucial point - energy and shape of the vortex solutions 
(pT|) - (jUnj) and (jUHjl - lfTOj) of the electric and magnet fields. To prove the finite- 
ness of energy of the vortex solutions (j^ - IU^ and (jUHjl - lfTOj) . we start with 
the equations for averaged in time balance of energy density of electrical and 
magnet field PH] : 



dW 



IGtt 



> dD* dD ^ 




(75) 



where D = Pu^ + AnPnUn is a sum of the linear induction and the nonlinear 
induction of the electrical field. In more of the books the calculations of the 
averaged energy of the optical waves in a dispersive media are worked out 
bearing in mind only the first order of slowly varying amplitude approxima- 
tion of the linear electrical induction. That is because the result comes to 
the old result of Brillouin 1921 for energy density of electrical field: 



lin 




(76) 



It is straightforward to show using this approximation that after integrating 
in space our vortex solutions they admit infinite energy. As seen from ()76|) . 
this approximation do not include the dispersion parameter /c", which is 
present in the envelope equations (j^ . (jUUI) as coefficient in front of the 
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second derivative in z direction. So there is no possibility to estimate the 
influence of k" in the energy integral. Our vector amplitude equations are 
obtained, using the second order approximation of Pim- The right expression 
of the energy density in this case requires to expand the linear electrical 
inductions in the energy integral also up to the second order, and to include 
this dispersion parameter in the energy integral. Using truncated form for 
first derivative in time of the linear induction of electrical field we have: 



dPiin . , . 7 cdA ick" d'^A 

^ = — o^(-o)'4 + -J^ + — g^. (77) 

From (ffSj) . (f77j) and complex conjugate of (f77|) we obtain: 



dW, 



lin 



dt 



c I ^dA* ^^dA 
V 'dt 



ick" d / -^^dA -^dA* 

J I A* A 

2 dt\ dt dt 



or 



Wiin = - 
V 



A 



+ 



ick" ( -^^dA -^dA* 



(79) 



where A* denote complex conjugate in time amplitude function. Rewriting 
the vortex solutions of electrical (jUT|) - (p3H|) and magnet (|nS|) - (f7n|) fields in not 
normed (dimension) coordinates and substituting in (f?^ . for the case of 
negative dispersion finally we obtain the next expression of the linear part of 
averaged energy density: 



Wun= [l + --c \ k" I Aujj A' , (80) 

where Aa; denote the spectral bandwidth of the vortices. We obtain one 
unexpected result: with the increasing of the spectral bandwidth of our solu- 
tions the linear part of energy density will decrease in the negative dispersion 
region. Conditions when the linear part of the energy density of electrical 
field is zero determine one critical spectral bandwidth Auc- 
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Aujc 



V + c 1 



(81) 



vc I k" 



Near to plasma frequency or some of the electronic resonances the dispersion 
parameter increases and values about: k" ~ 10~^^ — io~25cefc__ rjj^g critical 

^ cm 

spectral bandwidth in this case becomes: ~ 10^^ — 10^^ Hz. The envelope 
approximation request, that: ujq > AuJc- This condition shows that the linear 
part of energy density of vortices is zero only when their main frequency is 
situated in optical region, or regions, which frequency is greater than the 
optical and admit spectral bandwidth equal to Auc- The nonlinear part 
of averaged energy density is expressed in [22] and for the vortex solutions 
dnHl-dnni) it becomes: 



These results give the conditions for the finiteness of energy of the vortices: 
the spectral bandwidth of vortices to be equal to Auc- Integrating Wniin on 
the 3D space we obtained a finite value proportional to the main frequency 
ujq. The intensity part in this expression for the energy ()82|) is limited by the 
experimental condition (f7S|) and it is a constant. In this way we obtain the 
result, that the localized energy of the solutions increasing linearly with the 
increasing of the main frequency. The numerical investigation of stability 
for such type vortices, using split-step Fourier method, is provided in jS]. 
Comparing the numerical calculations for vortices with these, of a standard 
Gausian pulse, it can be clearly seen that the vortices propagate without any 
changing of their shape, whereas the Gaussian pulse, with the same initial 
amplitude, self-focuses rapidly in a much shorter distance. 

In additional, we also investigate the shape and the behavior of electrical 
and magnet field in the origin and infinity. To obtain the correct result, the 
real part solutions of the amphtude equations (jU^ - fjUH)) and (fUHj) - (jTOJ must 
be rewritten in the independent Cartesian coordinates x, y, z: 



-> ^ -> 2 (J II') -> J. -> I 

Wniin = n2A' A' +u:o^.A' A' . 



(82) 



A 




(83) 



x 



2;2 _j_ ^2 _j_ ^2 



A' 




(84) 



y 



2;2 _|_ ^2 _|_ ^2 
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The limits of these functions in origin and in infinity are investigated, using 
the Haine criteria for hmit of a muhidimensional function. It is straightfor- 
ward to show that the solutions (j83p - (j85p are odd functions, going to zero in 
infinity and admit very little value near to the origin. In the origin, using 
the same criteria, we find that the above functions have no limit. But as 
we must point out, the envelope function (jH^ - (jH3|) is only imagine functions 
around the real electric field, expressed in our case by the relation (|TTH) : 



Using the same criteria as in the previous case, it is seen, that the electrical 
field admit exact limit zero in the origin and infinity: 

lim^,y,z^o{Ei) = 
lim^^y^z^±oo{Ei) = 0- i = x,y,z. (89) 



6 Conclusion 



= A'^ sin {kQZ — UQt) 
Ey = A'y sin {koz — ujot) 
Ez = A'^ sin (k^z) — Uot. 



(86) 
(87) 
(88) 



The applied method of slowly varying amplitudes of the electrical and mag- 
net vector fields give us the possibility to reduce the nonlinear vector integro- 
differential equations to vector nonlinear differential equations of amplitudes. 
Here is the place to explain more some of the differences between the so- 
lutions, obtained by separation of the variables of the usual linear scalar 
Schredinger equation with potential depending only on 'r' (for example hy- 
drogen atom) and solutions of the vector version of NSE. For the linear 
Schredinger equations with potential, to higher order spherical functions 
correspond more higher order of radial spherical Bessel functions. While 
for vector version of NSE to higher order of spherical functions (£ = 1, 2, ..) 
correspond higher number of the fields components and higher value of local- 
ized energy. We have for all radial solutions the zero spherical Bessel function 
2^2-21 in nonlinear case. This is the reason to start with such special kind of 
complexification of our electrical field (jiup . It is important to highlight also 
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the equivalence in high frequency region between the hnear dielectric suscep- 
tibility of a cold plasma and this of a dielectric media. A plasma frequency 
of dielectric media in this region can determine also. The expression is equal 
to this of cold plasma wish precise a (dielectric) constant. This discussion 
shows two possible regions for observing optical vortices: 1. High frequency 
(transparency) region of cubic dielectrics or cold plasma. 2. Near to some of 
the electronics resonances of x^^^ media. Above mentioned investigations are 
provided only for the real part of linear and nonlinear susceptibility. That is 
why the first possibility, high frequency region, is more attractive for observ- 
ing of optical vortices. The complex part of linear susceptibility is significant 
near to electronic resonances and should be get in mind in a more detailed 
analysis. 
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